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A model for particle production at the deconfining phase 
transition in QCD is developed, as the semi-classical decay of 
a condensate for the Polyakov loop. In such a model, gener- 
ically particle production, as measured on an event-by-event 
basis, exhibits significant deviations from statistical behavior. 



The deconfmed phase of QCD at high temperature 
may be produced in the collisions of large nuclei at high 
energies, such as at RHIC and LHC. The usual picture 
of the deconfmed phase is of a weakly interacting gas of 
quasi-particles. While this picture applies (well) above 
the critical temperature, T c , it may not be very useful 
near T c . Recently, an effective theory which is better 
suited near T c was proposed; in a mean field approxima- 
tion, the free energy is due to a potential for the Polyakov 
loop § . 

In this Letter we propose that if the mean field theory 
is correct, particle production at the phase transition is 
driven by coherent oscillations of the Polyakov loop con- 
densate. This is very similar to particle production from 
preheating in inflationary models of the early universe . 
For both inflation and the mean field theory of [Q, the 
potential energy of the Polyakov loop dominates the total 
energy density, so that particle production can be esti- 
mated by semi-classical means. Particle production from 
a coherent source is observable by measuring fluctuations 
on an event-by-event basis, and differs from particle pro- 
duction by an incoherent, or statistical, source. 

Many of the features of our model are similar to 
particle production in Disoriented Chiral Condensates 
(DCC's) [HH]. Because pions are light, the chiral con- 
densate in a DCC has a small energy density relative to 
the total energy density. Thus the number of produced 
particles from the DCC is small relative to the total, and 
the particles are produced in a narrow region of phase 
space, on top of a large, incoherent background. In con- 
trast, in the present model essentially all of the energy 
of the deconfmed phase is going into oscillations of the 
Polyakov loop, and thereby into pions. We stress that 
while the Polyakov loop is treated classically, the pro- 
duction of pions is quantum-mechanical. 

We also contrast this with the conjectured critical end- 
point of the chiral phase transition in the plane of tem- 



perature, T, and (baryon) chemical potential, fi H[|. 
Through critical fluctuations in the sigma field, a crit- 
ical end-point will also generate non-statistical fluctua- 
tions. The critical end-point, however, only exists for a 
specific value of T and fj,; it is necessary to "tune" the 
collision parameters to reach this exceptional point. If 
the effective theory of jij applies to QCD, however, it 
should apply, about the critical temperature, for a large 
range of chemical potentials near zero; there is no need 
for tuning. Further, critical fluctuations J|,[| are neces- 
sarily largest about zero momentum. In contrast, in the 
present model the produced particles typically have mo- 
menta of order several hundreds of MeV, with significant 
fluctuations about that large momentum scale. Unlike a 
critical end-point, one should not cut on particles with 
low momentum; deviations from statistical behavior ap- 
pear in the average pion momentum. 

We begin by reviewing the mean field theory for the 
Polyakov loop Q, and then present an illustrative cal- 
culation of particle production. This calculation is not 
meant to be definitive, but merely to sketch how a clas- 
sical condensate for the Polyakov loop might decay. 

To understand the mean field theory, consider first an 
SU(N) gauge theory not for three colors, but in the limit 
of a large number of colors, N — ► oo. As first observed 
by Thorn [Q, in this case the free energy itself is an or- 
der parameter. For temperatures below the deconfining 
phase transition, confinement occurs, and the free energy 
is due exclusively to hadrons, such as glueballs. As color 
singlets, glueballs only contribute of order one to the free 
energy. We assume the usual, "quarkless" large N limit, 
in which the number of light quark flavors, Nf, is held 
fixed as N — > oo; then mesons also contribute of order 
one to the free energy. 

In contrast, above T c deconfinement occurs, and gluons 
contribute ~ N 2 to the free energy; the quark contribu- 
tion is suppressed, ~ N 2 (Nf /N). At very high temper- 
ature, by asymptotic freedom the behavior of the free 
energy can be computed perturbatively. 

There is a puzzle, however. The free energy is a gauge 
invariant quantity, and so at any temperature, should be 
describable exclusively in terms of gauge invariant exci- 
tations. While the glueballs change their nature with 
temperature, they remain the dominant gauge invariant 
excitations. But if color singlet glueballs can only con- 
tribute of order one to the free energy, what is the term 
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~ N 2 in the free energy due to? 

The only quantity which can provide such a contribu- 
tion is an expectation value for the Polyakov loop: 



£{x) = — tr | V exp | ig 



r l/T 



Ao(x,T)dT) ■ (1) 



tr is the matrix trace, V is path ordering, and g is the 
gauge coupling constant; Aq(x, t) is the time component 
of the vector potential in the fundamental representa- 
tion, at a spatial position x and euclidean time r. At a 
temperature T, the Wilson line is defined in imaginary 
time; we assume that it can be generalized to real time, 
although its detailed form is as yet unknown Q] . The full 
thermal Wilson line is an SU (N) matrix; its expectation 
value includes both color adjoint degrees of freedom, £ a , 
and a color singlet, I The £ a 's transform under the 
SU(N) gauge group, while the £'s only transform under 
a global Z{N) symmetry. The SU(N) Wilson lines, ~£ a , 
may be important in driving the transition first order 
when N > 3, but in a mean field approximation for the 
free energy, only the Polyakov loop, £, matter ||. (Fol- 
lowing standard convention, but contrary to p|, hence- 
forth we refer to the ts as Z(3) Polyakov loops, and to 
the 4's as SU(3) Wilson lines.) 

As the Polyakov loop is a color singlet, its expectation 
value, £q — (£), is gauge invariant. This vanishes in the 
confined phase, up to corrections ~ Nf/N ~ 1/N, and 
is nonzero above T c . At infinite temperature, by asymp- 
totic freedom we can ignore fluctuations in the gauge 
field, and £q — ► 1. 

Thus the term ~ N 2 in the free energy is due exclu- 
sively to the potential for the Polyakov loop condensate, 
V(£). The free energy of glueballs and other gauge invari- 
ant excitations contribute ~ 1, while quarks contribute 
- N 2 (Nf/N) ~ N. 

Of course N is not infinite in QCD, but three. At least 
for the free energy, perhaps large TV is a good guide to the 
behavior for N = 3. If so, then in the high temperature 
phase of QCD, the free energy is — in a precise sense 
- dominated by gluons, through the potential for the 
Polyakov loop. This is supported by Lattice data, which 
finds that the free energy is much smaller at low, than at 
high, temperature 111. 

There is one aspect of the deconfining phase transition 
for N = 3 which is very different from N = oo. General 
arguments and Lattice simulations || suggest that the 
deconfining transition is of first order when N > 4 (it may 
be strongly first order, although definitive results in the 
continuum limit arc lacking). In contrast, the deconfining 
phase transition is of second order for two colors j^] , and 
nearly second order for three colors jll| . Thus at least as 
far as the order of the deconfining transition in the pure 
glue theory is concerned, N — 3 is "near" N = 2, not 
N = oo. 



For N = 3, the order of the transition does change 
with the addition of dynamical quarks; we assume, how- 
ever, that even if the transition becomes crossover, that 
it remains nearly second order ffl. That is, we assume 
that Polyakov loops become light for some range of tem- 
peratures, which we define as T c . This is crucial to our 
analysis, because when £ becomes light, it oscillates with 
large amplitude about its minimum, and drives particle 
production. 

(As an aside, we note that the free energy need not 
be small in the low temperature phase of a gauge theory. 
Consider, for example, a generalized large N limit, in 
which Nf/N is held fixed as Nf and N — > oo. This 
large N limit is "quarky", as there are ~ N 2 color 
singlet hadrons. These hadrons generate a free energy 
which is of order N 2 ~ NNf ~ N 2 at any temperature. 
The Polyakov loop is also nonzero at any temperature, 
£q ~ Nf/N 7^ 0. Even though £ 7^ 0, perhaps the poten- 
tial for the Polyakov loop continues to dominate the free 
energy. However, in this case, even near T c the Polyakov 
loop is never light.) 

For three colors, the Polyakov loop £ is complex valued, 
and we take the potential 



V(£) = 



-|Ki 2 -|(^ 3 + (n 3 ) + J( 



|2\2 



64 T 



(2) 



£o(T) is the global minimum of V(£) at a temperature T. 
Again, we normalize £q — > 1 as T — > 00. 

An important feature is that because the Polyakov loop 
is the trace of a phase factor, it is a dimensionless field; 
the only scale to make up the correct powers of dimen- 
sion is the temperature. This accounts for the overall T 4 
in V(^). The coefficients 62, 63, and 64 are then fitted to 
agree with Lattice data for the pressure in the pure glue 
theory at T > T c ; the pressure for T < T c is taken to 
vanish. Of course the pressure in the pure glue theory is 
not the same as in QCD; for instance, the ideal gas value 
of the pressure changes with the addition of quarks. Lat- 
tice simulations have also shown that T c changes as well; 
with 2 + 1 flavors one finds T c - 180 MeV Lattice 
simulations have demonstrated, however, that even with 
dynamical quarks, the pressure, divided by the ideal gas 
value, is a nearly universal function of T/T c 0,0. We 
use this remarkable property in our fit, taking b2(T/T c ), 
and &3, from the pure glue theory The overall constant 
64 is rescaled by the ratio of the ideal gas terms between 
QCD, with three massless flavors, and the pure glue the- 
ory. 

The pressure in the pure glue theory for T > T c 
|Q is described with the constant values 63 = 2 and 
64 w 0.6061. In QCD, we take the same value of 63, 
and rescale 64 by the appropriate ratio of ideal gas val- 
ues, 64 w 0.6061 • 47.5/16. In the spirit of mean field 
theory, we take the same values for 63 and 64 below T c . 



2 



Below T c , an approximate expression for the string ten- 
sion is er(T) = 1.21^2 - 0.99T 2 /(0.41) 2 @, where a 
is the zero temperature string tension, = 1 GeV/fm. 
The mass of the Polyakov loop is then used to fix 62, 
m 2 = (<r(T)/T) 2 cx ~b 2 b 4 T 2 . 




FIG. 1. V(£)/T c 4 for temperatures T/T c 
to bottom), respectively. 



: 0.98, 1, 1.02 (top 



In Fig. [j] we show V(£) in the vicinity of T c , which is 
defined as the temperature were the two local minima of 
V(£) are degenerate. Notice that the potential changes 
extremely rapidly with temperature about T c . This is 
due primarily to the rapid change of the coefficient 62 
with temperature. This can be seen either above or be- 
low T c . Above T c , the "bump" in the trace of the energy 
momentum tensor, times 1/T 4 , is due to the variation 
of 62 it can also be seen from the extremely rapid 
decrease of the Debye screening mass, ~ 62, from 2T C to 
T c , by a factor of ten Below T c , 62 changes rapidly 
because of the quick rise in the string tension as the tem- 
perature decreases [Q. Physically, this occurs because 
the Polyakov loop is light at T c , and very heavy at zero 
temperature, me(0) > 1.5 GeV. 

Except near the critical temperature, the Polyakov 
loop is heavy, so its fluctuations, and particle produc- 
tion, are small. Only near T c does the Polyakov loop 
become light and generate pions. 

At T c , because there is a first order transition, there 
are two degenerate minima, at £ = and £ w 0.56. As 
can be seen from the figure, however, the barrier between 
these two minima is small. This is because the pure glue 
theory has a weakly first order transition when N = 3. It 
is so weak that this plays no part in our analysis. Thus 
whether the theory with dynamical quarks has a first 
order transition, or merely a crossover, is of little conse- 
quence; all that matters is that £ becomes light. 

(The rapid variation of the potential with tempera- 
ture, and the shallowness of the barrier between the two 
minima at T c , is in contrast to models of the chiral transi- 
tion Q| . For the chiral transition, the potentials change 
much more slowly with temperature, and the barriers to 
nucleation are significant; see, e.g., fig. (1) of 



The complete effective theory is described by adding a 
kinetic term for £, and coupling it to a chirally symmetric 
field, 6: 
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C = C^+'—\d^£\ 2 T 2 



V{£)--ct> 2 \£\ 2 T 2 . (3) 



(/> is the usual field describing pions and the sigma me- 
son, with 0(4) global symmetry. The model can also be 
extended to three flavors, with kaons, etc. 

Because the Polyakov loop is an effective field, normal- 
ized so that £0 — > 1 at high temperature, the coefficient 
of the kinetic term is not automatic. In the following 
analysis we require the coefficient of the kinetic term for 
the Polyakov loop, £, as it varys in time. We do not 
know this, so in (|J) we take the coefficient of the kinetic 
term for SU(3) Wilson lines, £ a , which vary in space [Q; 
we then assume a Lorentz invariant form to obtain the 
coefficient for time varying fields. While this is just a 
guess, taking a = g 2 /(4:ir) = 0.3 at T c gives a coefficient 
N/g 2 ~ 0.8, which is about unity and so reasonable. 

The coupling constant for the Polyakov loop with 
mesons, h 2 , could be fixed by knowing £q(T) and the 
meson masses above T c . Since £q(T) is not known [Q, 
our estimate is imprecise. Comparing to Lattice data of 



Gavai and Gupta for m w (T) at T > T c |15|, m^/T = 2tt 
at T — > 00 gives h 2 « 29; at T/T c = 2.9, m^/T = 4.29, 
and the Hartree relation rr£ = A 2 (T 2 /2 - v 2 ) + h 2 T 2 \£\ 2 
(which follows from eqs. [|||), yields h 2 = 9. In the actual 
calculations described below we employed h 2 = 22 as an 
intermediate value. 

The lagrangian for the <j> field is standard: 



b 2 -v 2 ) 2 



Ha 



(4) 



At T = chiral symmetry is spontaneously broken, and 
the expectation values are (a) = f w , (7?) = 0. Further, 
the PCAC relation gives H = f^m 2 . , were = 138 MeV 
is the pion mass in the vacuum state at T = 0. Then 
v' 2 = f 2 — m 2 /A 2 , and assuming a coupling A 2 = 20 
yields a sigma mass in vacuum of m 2 — 2A 2 / 2 + m 2 w 
(600 MeV) 2 . 

There are many other terms which we could have in- 
cluded in the effective lagrangian. In particular, we really 
should have included terms linear in the Polyakov loop. 
In the pure glue theory, such terms are prohibited by 
the global Z(3) symmetry of 't Hooft, but it is expected 
to occur with dynamical quarks ]l6[| , A term linear in 
the Polyakov loop will change the first order transition 
(in the pure glue theory) to crossover. We have ignored 
terms linear in the Polyakov loop because the free energy 
is small below T c , so they must be small. 

In this vein, we note the results of Digal, Laermann, 
and Satz jL7|. They obtained Lattice results for two 
flavors of light, dynamical quarks. As expected for a 
second order chiral phase transition, about T c there is 



3 



a peak in the chiral susceptibility, which becomes more 
pronounced as the quark mass decreases. For the range 
of quark masses studied, however, the peak in the sus- 
ceptibility for the Polyakov loop also continues to grow. 
A term linear in the Polyakov loop should eventually cut 
off the divergence in the susceptibility for the Polyakov 
loop; thus their results also suggest that in the effective 
lagrangian, the effects of all terms linear in £ (times one, 
</> 2 , etc.) are small. 

To compute particle production we take a simple dy- 
namical picture. We assume an instantaneous quench 
from an initial temperature T + > T c , to a temperature 
T_ < T c H . We stress that because the ^-potential varies 
so rapidly with temperature around T c , the assumption 
of a quench appears natural. We then evolve the field, 
given its value for the minimum at T+, with the potential 
at T_ . The £- field is no longer at a minimum at T_ , and 
so it rolls down the potential, and oscillates around the 
new minimum at £q = 0. These oscillations in £ couple 
to 4>, and thereby produce pions and sigmas. For the 
purposes of illustration we only include the production 
of pions, ignoring that of sigma mesons. 

We perform the standard decomposition of the pion 
field in terms of creation and annihilation operators times 
adiabatic mode functions n k (t) 0Jll|. The equation of 
motion for the mode functions is 



dt 2 



= -n 2 k w k (t) , 



n 2 =k 2 +m 2 n + \ 2 (ir 2 (t))+h 2 \£\ 2 T 2 



(5) 
(0) 



The effective pion mass is calculated assuming Hartree 
type factorization, 



7T 2 (t,x)^(n 2 (t)) 
which requires a subtraction at t = Jl£ 
d 3 k _ , 



(7) 



(*'(*)) 



(27t) 3 

d 3 k f N k {t) 
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(2tt)3 \Q k (t) 2fi k (t) 2n k (t = Q) 



(8) 



d-n = 3 counts the number of isospin states. Inserting 
a cutoff A, which is assumed to be larger than all other 
mass scales, gives 



d w 



Co log , , 



4A 



c 2 log 







4A : 
C 
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d 3 k N k 
(2ir) 3 n k 



(9) 



where c 2 , = to 2 + h 2 \£(t = 0)| 2 T 2 and c 2 = m 2 + 
A 2 (7? 2 (t)) + h 2 \£(t)\ 2 T 2 . In the above equations, N k de- 
notes the occupation number of a mode with given isospin 
and with momentum k, subtracted for vacuum fluctua- 
tions, 



A4 = 



(10) 



As initial conditions we consider the case that the state at 
t = contains vacuum fluctuations only, such that -K k (t — 

0) = y/l/2Sl k (t = 0), dw k (t = 0)/dt = -iy/Q k (t = 0)/2, 

and N k (t = 0) =0. 

Similarly, in the classical equation of motion for the 
Polyakov loop, we include the backreaction from the 
produced pions in the Hartree approximation, replacing 
7r 2 (i, x) by (n 2 (t)) as given in (||). We neglect fluctuations 
of £ here, which should be a reasonable approximation as 
long as the energy of the produced pions is well below 
that for the Polyakov loop at t = 0. Of course, a more 
refined treatment is necessary in order to trace the time 
evolution to the point where the coherent oscillation of £ 
has dissipated its entire energy into pions, since produced 
mesons will scatter off the Polyakov loop condensate and 
cause it to decohere. 

We performed numerical simulations for the following 
values of the parameters: A = 10/^ ~ 1 GeV, T + — 1.02, 
T_ = 0.975. In preheating after inflation the coupling 
constant is very small, ~ 10~ 12 , and the inflaton field 
oscillates many times. Particle production occurs con- 
tinuously during these oscillations, and resonance bands 
develop for modes in phase with the driving inflaton field. 
In the present model, both the coupling of the Polyakov 
loop to pions, h 2 = 22, as well as the self coupling of the 
pions, A 2 = 20, are large. Consequently pion production 
happens rather quickly, within a few oscillations. 




k/f. 



FIG. 2. Occupation numbers of the pion modes at times 
t = 4n/3/,r, for n = 1..3. 

In fig. ^ we show the logarithm of the produced pion 
occupation number versus momentum. The distribution 
exhibits less pronounced resonance bands than in the 
case where the contribution of fluctutations to £l k is ne- 
glected (compare also fig. 4 and fig. 5 in ]l9| ). That is 
because X 2 (n 2 (t)) is time dependent, and shifts the reso- 
nance bands to other k modes. 
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FIG. 3. Pion density as a function of time. 



Integrating over momentum, 



Pit 



2tt 2 



dk k 2 M k 



(11) 



we obtain the total density of produced pions. As shown 
in fig. ||, the pion density increases very rapidly: Para- 
metric resonance lets the pion density increase exponen- 
tially over the depicted time interval. 

The total pion density at time t = 4//^ after the 
quench is approximately p v = 0.35/fm 3 . From p„, the 
total number of pions per unit rapidity, dN/dy, is esti- 
mated as follows. Longitudinal length in the beam direc- 
tion is a scale factor, a, times length in rapidity. For one 
dimensional Bjorken expansion P0| , a = t, where r is the 
proper time. If the expansion is isentropic, and s is the 
entropy density, then sa is a constant; taking s ~ T 3 , 
a rough estimate of a(T c ) at RHIC is ~ 10 fm. Then 
dN/dy = p w 7rii 2 a(r c ). For a nucleus with A ~ 200, 
even without transverse expansion R ~ 7 fm, so that 
dN/dy « 1300/9^ fm 3 sa 460. If transverse expansion 
increases R by 50%, dN/dy doubles and agrees roughly 
with data for central collisions at BNL-RHIC |pl| . 




FIG. 4. Mean energy per pion (in a given event and do- 
main) as a function of time. 

Fig. |5] shows the average energy per produced pion as 
a function of time. Unlike a thermal bath, the energy of 



all pions oscillates in phase with the background £-field, 
through the terms (n 2 ) and h 2 \£\ 2 . Because there is a 
lot of energy in the £-field, the pion energy is relatively 
large, and there are large fluctuations about the aver- 
age. From Fig. |1| the time-averaged root-mean-square 
(RMS) fluctuation of the average pion momentum is 
y/(p 2 ) — (p) 2 /(p) ~ 10.3%. We have taken the ensem- 
ble average, i.e. the average over events, to be given by 
the time average from t = up to t = 4//^. 

This is our principal result. In a single domain, as for 
figs. U and |], decay of a classical Polyakov loop conden- 
sate generates large fluctuations. In this picture, pions 
are produced in a pulse near T c , from the "ringing" of 
the Wilson line into pions. 

There are two effects which act to reduce the effect 
dramatically. First of all, experimentally it is neces- 
sary to average over many domains. The size of each 
domain is determined by how the quench ends, and in- 
volves both quantum mechanical and semi-classical pro- 
cesses. A lower bound for the size of a single domain is 
presumably given by the correlation length of the ^-field 
at T c , £ ~ l/mg(T c ). From our estimates above, because 
I becomes light at T c , this correlation length is relatively 
large, £ ~ 1 fm. Assuming one samples experimentally 
the full transverse size of the nucleus, and over one unit of 
rapidity, gives N D = 7ri? 2 a(T c )/(47r£ 3 /3) ~ 300 domains. 
While perhaps semi-classical methods are suspect when 
each domain gives ~ 1.5 pions, if we average over 300 
domains, the RMS fluctuation of the average momentum 
per pion is reduced to ~ 10.3% /^/Njj = 0.6%. Averaging 
over many domains also acts to smooth out the momen- 
tum distribution of fig. to leave a distribution which is 
closer to exponential. 

While fluctuations ~ 0.6% are small, they are experi- 
mentally measurable. For central collisions at the CERN- 
SPS, y/s = 18^4 GeV, the NA49 experiment compared 
the following two distributions of the average pion mo- 
mentum ]22||B3||. The first was the event-by-event dis- 
tribution of the actual data, the variance of which reflects 
intrinsic physical correlations in the event. The second 
is the same distribution in mixed events, constructed by 
picking particles at random from different events. The 
variance of mixed events is determined by purely sta- 
tistical fluctuations, such as those due to finite particle 
number. To within 0.05%, however, the RMS fluctua- 
tions of the two distributions are the same. 

Apparently, our model does not apply at SPS ener- 
gies. One possibility is that the model is wrong. An- 
other is that too large a region in rapidity, etc. was 
averaged over. It would be interesting to bin the data 
in increasingly small bins in rapidity, until one is limited 
by the usual statistical uncertainty, ~ 1/y N& y , where 
NAy — (dN/dy) Ay is the number of particles in a rapid- 
ity bin of width Ay. 

It is also possible that fluctuations at the SPS are 
washed out by scattering in a hadronic phase. We as- 
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sume that particles produced at T c essentially flow to 
the detectors without further collisions, which is of course 
an idealized scenario. Collisions in a hadronic phase act 
to erase intrinsic fluctuations generated at T c , and leave 
purely statistical fluctuations. We leave a detailed in- 
vestigation of this point for future work. It is not clear, 
however, how easy it is for hadronic scattering to wash 
out fluctuations at T c , due to expansion inherent in a 
heavy ion collision. Pions are nearly Goldstone bosons, 
whose scattering is suppressed at low energy by powers 
of momentum. For example, in model calculations |l9| ], 
even though the self coupling of pions is large, pions can 
scatter into the zero mode of the pion field more rapidly 
than they can scatter out. 

Our estimate of dN/dy is just to tune the parame- 
ters of the model, but fluctuations in mean momentum 
~ 0.6% should be a qualitatively correct prediction. In- 
deed, one can easily imagine how the fluctuations could 
even be larger. If the average size of the domains, £, is 
not ~ I but ~ 2 fm, then in one unit of rapidity there are 
only 40 domains instead of 300. Each domain generates 
about ten pions, and the RMS fluctuations of the average 
momentum per pion increases to ~ 10.3%/a/40 = 1.7%. 
Notice, too, that as particle production increases expo- 
nentially in time, the distribution is dominated by those 
domains which last for the longest time. 

If indeed hadron production at the confinement tran- 
sition is dominated by the decay of a Polyakov loop con- 
densate, we also expect nonstatistical (intrinsic) fluctu- 
ations of the pion multiplicity from event to event. For 
example, the RMS fluctuation of the multiplicity during 
the last 1 fm/c depicted in fig. [| is ~ 24% in a single 
domain; for 300 domains, this gives fluctuations ~ 1 .4%. 

Recently, fluctuations of the electric charge were ana- 
lyzed Based upon a quasi-particle picture, it is ar- 
gued that that RMS fluctuations of the charge, measured 
in a rapidity window Ay ~ 1, should be much smaller in 
a quark-gluon plasma than in a pion gas. Thus, if a 
thermal and locally neutral plasma were produced, and 
the pions emerged from the hadronization of that plasma 
with small relative rapidity, the charge fluctuations in the 
final state remain small: given a charged pion in Ay, the 
probability for its partner with opposite charge to have 
rapidity within Ay is large. In the present model, on 
the other hand, pions produced from oscillations of the 
electrically neutral Polyakov loop are spread over a ra- 
pidity interval of ~ ±2 (see fig. ^|), and the number of 
7r + and 7r~ in a given rapidity window Ay ~ I fluctuates 
independently. Thus, the RMS fluctuations of the charge 
should be large, i.e. like that of hadrons, not quarks. 

The fluctuations are dominated by the behavior of the 
system as it goes through T c . To first approximation, it 
does not matter how hot the system is initially. Con- 
sequently, the magnitude of the fluctuations should be 
approximately similar in collisions at RHIC and at the 
LHC. The only change with increasing the energy of the 



colliding nuclei is a change in the scale factor a(T c ), which 
produces a higher multiplicity at higher energy. As a 
constant number of pions is produced per domain, the 
width in rapidity should be chosen so that there is the 
same number of particles per rapidity bin. Higher ener- 
gies can help, albeit indirectly: transverse expansion may 
be larger at higher energies, which prevents rescattering 
in the hadronic phase from washing out fluctuations in- 
duced at T c . 

We have made numerous approximations in the present 
analysis. A more careful study is certainly possible. The 
potential for the Polyakov loop in QCD, and the coupling 
to pions, can be computed from numerical simulations 
on the lattice, both without and with dynamical quarks. 
Many other improvements can also be made: introduc- 
ing realistic models for the space-time dependence of the 
temperature, the effect of expansion in the equations of 
motion for t and <f>, etc. 

We have not attempted a better analysis here because 
STAR and other detectors will soon decide if heavy ion 
collisions at RHIC do or do not display intrinsic fluctua- 
tions on top of purely statistical fluctuations. We eagerly 
await this analysis. 

Note added: After this paper was submitted for pub- 
lication, many results were announced at Quark Matter 
2001. CERES reported data from the SPS indicating 
fluctuations in the average transverse pion momentum 
- 3% ||. At RHIC, STAR found that the same fluc- 
tuations appear to be much larger, ~ 8% j26|. STAR 
also reported that pion interferometry appears to indi- 
cate "explosive" behavior |^6|,^7j . We have recently sug- 
gested how the Polyakov loop model might generate ex- 
plosive particle production and large fluctuations pjl . 



ACKNOWLEDGMENTS 

A.D. acknowledges support from a DOE Research 
Grant, Contract No. DE-FG-02-93ER-40764; R.D.P. 
from DOE grant DE-AC-02-98CH-10886. R.D.P. thanks 
the Heavy Ion Group at Wayne State University for dis- 
cussions, including R. Bellwied, S. Gavin, C. Pruneau, 
and S. Voloshin. 



[11 R. D. Pisar ski, Phys. Rev. D62, 111501 (2000); hep 



ph/0101168 



[2] J. H. Traschen and R. H. Brandenberger, Phys. 
Rev. D42, 2491 (1990); L. Kofman, A. Linde and 
A. A. Starobinsky, Phys. Rev. D56, 3258 (1997). 

[3] A. A. Anselm, Phys. Lett. B217, 169 (1989); AA. 
Anselm and M.G. Ryskin, Phys. Lett. B266, 482 (1991); 



6 



J.-P. Blaizot and A. Krzywicki, Phys. Rev. D46, 246 

(1992) ; J.D. Bjorken, Int. J. Mod. Phys. A7, 4189 (1992). 
K. Rajagopal and F. Wilczek, Nucl. Phys. B404, 577 

(1993) . 

M. A. Halasz, A. D. Jackson, R. E. Shrock, 
M. A. Stephanov and J. J. Verbaarschot, Phys. Rev. 
D58, 096007 (1998). 

M. Stephanov, K. Rajagopal, and E. Shuryak, Phys. Rev. 
Lett. 81, 4816 (1998); Phys. Re v. D60, 114028 (1999); 
B. Berdnikov and K. Rajagopal, |iep-ph/9912274 . 



C. B. Thorn, Phys. Lett. B99, 458 (1981); R. D. Pisarski, 
Phys. Rev. D29, 1222 (1984). 

S. Ohta an d M. Wingate, Nucl. Phys. Proc. Su ppl. B73, 



435 (1999); |hep-lat/9909125| ; [hep-lat/0006016 



J. Engels and T. Scheideler, Phys. Lett. B394, 147 

(1997); Nucl. Phys. B539, 557 (1999). 

F. Karsch, Nucl. Phys. Proc. Suppl. 83, 14 (2000). 

F. Karsch, E. Laermann and A. Peikert, Phys. Lett. 

B478, 447 (2000). 

B. Beinlich, F. Karsch, E. Laermann and A. Peikert, Eur. 
Phys. J. C6, 133 (1999). 

O. Kaczmarek, F. Karsch, E. Laermann and M. Lutge- 
meier, Phys. Rev. D62, 034021 (2000). 
O. Scavenius, A. D umitru, E. S. Fra ga, J. T. Lenaghan 
and A. D. Jackson, hep-ph/000917l|. 



R. V. Gavai and S. Gupta, |hep-lat/00040li 



T. Banks and A. Ukawa, Nucl. Phys. B 225, 145 (1983) 



S. Digal, E. Lae rmann, and H. Satz, hep-ph/0007175 



hep-lat/0010046 



D. Boyanovsky, H. J. de Vega, R. Holman and J. F. Sal- 
gado, Phys. Rev. D54, 7570 (1996); D. Boyanovsky, 
D. Cormier, H. J. de Vega, R. Holman, A. Singh and 
M. Srednicki, Phys. Rev. D56, 1939 (1997). 

A. Dumitru and O. Scavenius, Phys. Rev. D62, 076004 
(2000). 

J. D. Bjorken, Phys. Rev. D27, 140 (1983). 

B. B. Bac k et al. [PHOBOS Collaboration], |hcp" 



ex/0007036 



H. Appelshauser et al. [NA49 Collaboration], Phys. Lett. 
B459, 679 (19 99). 



H. Heiselberg, nucl-th/0003046, and references therein 



M. Asakawa, U. Heinz and B. Muller, Phys. Rev. Lett 
85, 2072 (2000); S. Jeon and V. Koch, Phys. Rev. Lett 



nucl- 



85, 2076 (2 000); H. Heiselberg and A. D. Jackson.. 
th/000602i| M. Bleicher, S. Jeon and V. Koch, frep- 
ph/000620 



H. Appelshauser, to appear in the proceedings of Quark 
Matter 2001. 

J. Harris, S. Panitkin, and J. Reid, to appear in the pro- 
ceedings of Quark Matter 2001. 

F. Laue and S. Panitkin, to appear in the proceedings of 
QM 2001. 



A. Dumitru and R. D. Pisarski, |hep-ph/0102020[ , to ap 
pear in the proceedings of QM 2001. 



7 



